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ABSTRACT 

We present an analysis of high-resolution hydrodynamical N-body simulations of cou- 
pled dark energy cosmologies which focusses on the statistical properties of the trans- 
mitted Lyman-a flux in the high-redshift intergalactic medium (IGM) . In these models 
the growth of the diffuse cosmic web differs from the standard ACDM case: the density 
distribution is skewed towards underdense regions and the matter power spectra are 
typically larger (in a scale dependent way). These differences are also appreciable in 
the Lyman-a flux and are larger than 5 % (10%) at z = 2 — 4 in the flux probabil- 
ity distribution function (pdf ) for high transmissivity regions and for values of the 
coupling parameter /3 = 0.08 (/3 = 0.2). The flux power spectrum is also affected at 
the ~ 2% (~ 5 - 10%) level for P = 0.08 (/3 = 0.2) in a redshift dependent way. We 
infer the behaviour of flux pdf and flux power for a reasonable range of couplings and 
present constraints using present high and low resolution data sets. We find an upper 
limit /3 < 0.15 (at 2a confidence level), which is obtained using only IGM data and is 
competitive with those inferred from other large scale structure probes. 
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1 INTRODUCTION 



, The observational evidence of the present accelerated cos- 
' mic expansion represents one of the major challenges to 
our understanding of the Universe. The standard ACDM 
■ cosmological model identifies the origin of this acceleration 
' with a cosmological constant term in the field equations of 
General Relativity. However, this interpretation suffers of 
extremely severe fine-tuning problems, and possible alter- 
native explanations of the accelerated expansion in terms 
of a Dark Energy (DE) dynamica l field have been pro- 
posed (as e.g. by IWet tcrich 198^: iRatra fc Peebled [19881 : 
lArmendariz-Picon et al.ii200Q ). Among these, particular at- 
tention has been recen t ly devoted to coupled DE models 
(cDE ) jWcttcrich 1 19951 : lAmendolal l2000l : iFarrar fc Peebles! 
I2OO4I : iBaldi 2010) where a direct interaction between the 
DE field and Cold Dark Matter (CDM) particles determines 
peculiar fe atures i n the background expansion of the Uni- 
verse {e.g. lAmend ola 20 00), in the evolution of linea r den- 
sity perturbations {e.g. iDi Porto fc Amendolal 120081 ') . and 
even in the nonli near dynamics of collapsed s tructures at 
small scales {e.g. iBaldi et al.ll2010l : lBaldill2010l '). It is there- 
fore of crucial importance in the present cosmological in- 



vestigation to devise observational tests capable to distin- 
guish between the standard ACDM cosmol ogy and alter- 
nativc DE models as e.g. the cDE scenario (jHonorez et al.l 
2010; Baldi & Pcttorino 2010|). In the present work, we ex- 
plore the possibility to use the observed properties of the 
diffuse baryonic matter at high redshifts as a direct probe 
to test and constrain cDE cosmologies. 

Standard cosmological models based on cold dark mat- 
ter plus a cosmological constant predict that most of the 
baryons at high redshift are in a diffuse form, the Inter- 
galactic Medium (IGM), and fill a significant portion of the 
Universe, giving rise to the so-called cosmic web; a network 
of median fiuctuated filaments interconnecting galaxies and 
tracing the underlying dark matter distribution. A great 
progress in the study of the IGM has been recently made 
thanks to the large data sets available and in particular 
high resolution quasar (QSO) spectra ( Ultra Violet Echelle 
Spectrograph or HIRES) and the low resolution Sloan Digi- 
tal Sky Survey (SDSS) QSO spectra: the present limitations 
appear to be of systematic nature rather than statistical. 
High and low-resolution Lyman-a QSO spectra of distant 
sources are thus very useful in characterizing the proper- 
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ties of the underlyin g mass density field at z = 2 — 6 alori] 



theline-of-sight (e.g. Isi fc Davidsenlll997l:[Croft et al| 200 



dp 

m 



IViel et al.ll2004l : iMcDonald et al.ll2005l : iMeiksinI 120091 ) and 
are now routinely analyzed to reconstruct the matter distri- 
bution in three-dimensions. The dynamical information on 
the growth of structures is however convolved non-linearly 
with other physical effects determining the shape of the 
Lyman-Q absorption lines in redshift space (thermal history, 
peculiar velocities, etc.), and also observational procedures 
(continuum fitting, removing of the metal lines and strong 
absorption systems, etc.) need to be prope rly modeled and 
marginalized over (see e.g. lMcDonaldll200(i ) . All these issues 
make it quite difficult to use the IGM as a cosmological tool 
since its structures have to be modeled with hydrodynami- 
cal simulations that incorporate the most relevant physical 
processes. Nevertheless, IGM Lyman-Q data in combination 
with other probes provide the tightest constraints to date 
on the neutrino mass fraction and on the coldness of CDM 
and suggest a higher r.m.s. value for the amplitude of the 
matter power than th at obtained from cosmic mi crowave 
background data [e.g. ISehakll2005l : IViel et al.ll2007^ . 

In this work we perform the first high-resolution hy- 
drodynamical simulations of cDE models with gas cooling 
and star formation in order to check whether one can use 
the IGM to detect the difi'erent growth of cosmic structures 
predicted by these models compared to the standard ACDM 
scenario. To do so, we make use of the modified versi on of the 
paral lel hydrodynamical N-body code GADGET-2 (jSpringell 
l200a ) specifically devised to include the effects of cDE mod- 
els and already presented in lBaldi et al.l (|201G| ). 

The layout of this Letter is as follows: in Sec. 2 we 
describe the cosmological models investigated and the hy- 
drodynamical simulations, in Sec. 3 we present the results 
in terms of 1-pt and 2-pt matter and flux statistics, in Sec. 4 
we compute the constraints on the DE-CDM coupling pa- 
rameter /3 using the SDSS flux power and the flux pdf of 
UVES spectra and in Sec. 5 we draw our conclusions. 



2 HYDRODYNAMICAL SIMULATIONS OF 
COUPLED DARK ENERGY MODELS 

We consider cDE cosmologies where the accelerated expan- 
sion of the Universe is driven by a DE scalar field <j) which 
interacts with the CDM fiuid by directly exchanging energy 
according to the equations: 



p'c + iT-Lpc 
P^ + ZHpij, 



-Pcj)' Pc 

+/3<j}'pc 



(1) 
(2) 



where a prime denotes a derivative w.r.t. conformal 
time and H is the conformal Hubble function. The con- 
stant parameter /3 fully specifies the interaction and 
determines the strength of the DE-CDM coupling^. 
The cDE models described by Eqs. Hll2p have been 
widely studied i n the li terature (see e.g. WetterichI 
1 19951 : lAmendolal l200d : iPettorino fc Baccigalupil l2008l : 



^ Note that we have used units in which Mpi = (SttG)"^''^ = 1 
and that the definition of the coupling /3 differs by a factor ^3/2 
from the one used in some of the literature. 



IWintergerst fc Pettorinoll2010l : iBaldi et aLlbOld and refer- 
ences therein) to which we refer the reader for an exten- 
sive discussion of the main features of these models. For 
the analysis carried out in this Letter, the effect of primary 
interest is the enhanced growth of CDM density perturba- 
tions arising in cDE cosmologies due both to the long-range 
fifth-force acting between CDM coupled particles (which 
attract each other with an effective gravitational constant 
G = Gjvll -I- 2/3^], where Gn is the standard Newtonian 
value) and to the additional velocity-dependent acceleration 
of coupled particles a„ oc pcj>v which follows from momentum 
conservation in cDE models. 

For our analysi s we rely on simul ations run with the 
modified version by I Baldi et al.l (|2010l l of the parallel hy- 
drodynamical N-body code GADGET-2 based on the con- 
servative 'entropy-fornmlation' of Smoothed Particle Hydro- 
dynamics (SPH) (i^Springcl 2005). They consist of a cosmo- 
logical volume with periodic boundary conditions filled with 
an equal number of dark matter and gas particles. Radia- 
tive cooling and heating processes were followed for a pri- 
mordial mix of hydrogen and helium. The star formation 
criterion simply converts in coUisionless stars all the gas 
particles whose temperature falls b elow 10^ K and whose 
density contrast is larger than 1000 (|Viel et al.l (|2004 ) have 
shown that the star formation criterion has a negligible im- 
pact on the fiux statistics considered here). More details on 
the gas cooling and on th e Ultra Violet background can be 
found in I Viel et al] (|2004l '). 

The cosmological reference model corresponds to a 
'fiducial' ACDM Universe with parameters, at z = 0, ilm ~ 
0.26, Qa = 0.74, fib = 0.044, n^ = 0.963, Ho = 72 
km s~^ Mpc~^ and as — 0.7 96, consistent with t he re- 
suhs of WMAP 5-years data JKomatsu et al.l bOOgl ). The 
cDE models considered in our analysis are two scalar field 
models with a Ratra-Peebles (Ratra fc Peebles 1988) self- 
interaction potential, U{<j>) — Uo4>~" , with a = 0.143 
and with couplings (5 = 0.08 and P — 0.2 , alrea dy stud- 
ied in iMaccio et al.l (|2004l ) and iBaldi et all l|2010t ). and la- 
belled "RP2" and "RP5", respectively. These acronyms re- 
fer to the type of self-interaction potential (where "RP" 
stands for the Ratra-Peebles potential) and to the strength 
of the coupling in a scale ranging from /3 — 0.04 (which 
would appear as RPl) to /? = 0.2 (our RP5 model). 
Other numerical simulations for some related models of 
interacting DE or modified gravity - but without hy- 
drody n amics - have been presen t ed by e.g. iNusser et al-l 
(l2005h: iKesden fc Kamionkowskil feOOCT: ' Schmidtl l|2009l )i 
iKeselman et al.l l|2010l ): lLi fc Barro w (2010). 

Initial conditions for all the runs are realized by displac- 
ing particles from a cartesian grid according to Zel'dovich 
approximation in order to obtain a particle distribution 
with the desired spectrum of density fluctuations. The ini- 
tial shape of the power spectrum for cDE models is slightly 
blue-tilted with r espect to the fiducial ACDM scenario (see 
iBaldi et al.ll2010l . for a description of the tilt in cDE mod- 
els) and this effect is taken into account in our runs. The 
amplitude of the initial {z — 60) density fluctuations in the 
different simulations is then rescaled with the appropriate 
growth factor in order to normalize all the linear power spec- 
tra to the same erg at 2 = 0. It is also important to stress 
that the cDE and ACDM simulations have been performed 
with the same random phases in the initial conditions and 
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Figure 1. Left: The IGM density pdf of tlie RP5 (tliick curves), and the ACDM (thin curves) models at 2 = 2.2,3,4.2 (red dashed, 
black continuous, blue dotted curves, respectively). Right: Matter power spectra of the ACDM and the RP5 models. 



the same set of cosmological parameters at ^ = 0: thereby 
the differences that we will highlight will be exclusively due 
to the different couplings and not to other parameters. We 
have used 2 x 400 dark matter and gas particles in a 60 ^"'^ 
comoving Mpc box for the flux power to sample the scales 
probed by the SDSS. The gravitational softening was set to 
5 h~^ kpc in comoving units and the mass per DM particle 
is 2 X IO^Mq/Zi. The snapshots analysed are in the redshift 
range z — 2.2 — 4.2 since these are the redshifts at which the 
flux power is measured by SDSS and 1000 quasar spectra are 
extracted for each redshift bin. Noise properties are added 
to the spectra in order to reproduce the observed data. 



3 RESULTS 

3.1 Impact on density and flux statistics 

We first consider the matter fields as extracted from the 
simulations. In Figure [l] (left panel) we show results for the 
IGM density probability distribution functions between cDE 
models and ACDM at 2; = 2.2, 3, 4. The red, blue and black 
curves refer to the three different redshifts {z = 2.2, 3, 4, 
respectively) while the thick ones represent the results for 
RP5 (/3 = 0.2) and the thin ones for ACDM (for clarity we 
do not plot the RP2 model in these figures). It is clear that 
it is more likely to encounter low density regions along the 
line-of-sight in a cDE model rather than in a standard cos- 
mology: i.e. the gas distribution of a cDE model is skewed 
towards regions that are less dense than the mean. While 
the differences at 2: = 4 for the RP2 model are about a fac- 
tor two for p / < p >= 0.1, they rise to a factor seven for 
the RP5 model. A visual inspection of some one-dimensional 
gas density fields confirms this trend and shows that in cDE 
models the underdense regions are usually emptier than in 
the corresponding ACDM case. There are also some differ- 
ences at large densities but we do not focus on those here 
since their volume filling factor is much smaller than less 
dense regions. 

In the right panel we plot the matter power spectra 



for the same two models shown in the left panel: the non- 
linear power spectrum of cDE models is tilted and usually 
larger than for the ACDM model. We recall that the power 
spectra are normalized to reproduce the same erg at 2: = 
and thereby their growth is different (see lBaldi et al.ll2010l ). 
The pivot scale is at about 0.25 h/Mpc, which corresponds 
to (roughly) a scale of about 0.0025 s/km in redshift space 
for the corresponding one dimensional flux power. At scales 
1-2 h/Mpc the matter power for RP5 (RP2) is about 25% 
(2%) larger than the ACDM matter power (we have again 
plotted here only the ACDM and the RP5 models for clarity 
reasons) . 

We note that the inferred linear matter power spec- 
tr um amplitude from the analysis of the SDSS flux power 
bv lMcDonald et al.l ((20051) is recovered with a la error bar 
of about 15%, thus we naively expect, even without per- 
forming any analysis in terms of flux power, that this k— 
dependent increase in the matter power can be constrained 
and the RP5 model ruled out from this data set alone. 

We now turn our attention to the mock sets of QSO 
spectra that have been extracted from each snapshot and 
normalized to reproduce the same (observed) effe ctive opti- 
cal de pth Tcs{z) = — In < F{z) > as estimated bv lKim et al.l 
120071) . The scaling factors that are used in order to normal- 
ize the optical depth differ by less than 4% between cDE 
and ACDM models. 

In Figure [2] we plot the 1-pt and 2-pt distribution func- 
tions in terms of flux pdf (left panel) and flux power (right 
panel). The color coding is the same as in the previous fig- 
ures. The differences that are present for the low-density tail 
of the gas pdf show up in regions of high transmissivity (flux 
values ~ 1). This demonstrates that cDE voids are emptier 
of gas and contain less neutral hydrogen than in a corre- 
sponding ACDM universe. Our findings therefore suggest 
that cDE models might alleviate the discrepancy between 
simulations and observations in voids, thereby pro viding a 
possible solution to the so called "Void Problem" (jPeebled 
I2OOII : iKeselman et al.l I2OI0I ) . It is also intrigu i ng th at from 
the analysis of SDSS data [McDonald et alj (|2005l ) find a 
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Figure 2. Left: Ratio between the flux pdf of the RP5 (thick curves), RP2 (thin curves) and the ACDM models at 2 = 2.2,3,4.2 (red 
dashed, black continuous, blue dotted curves, respectively). Right: Percentage differences between the flux power of the RP2 and RP5 
models and that of ACDM. The shaded areas in both panels represent the statistical error at 2 = 2.94 for the high resolution PDF data 
of Kim et al. (2007) and for the 2 = 3 bin of the SDSS QSO flux power as computed by McDonald et al. (2006). 



structure growth at z = 3 that is faster than that predicted 
by ACDM, although with a relatively low statistical signifi- 
cance: this could be better reproduced by cDE rather than 
standard cosmologies. 

We also note that there is an increase in regions of 
low transmissivity, although this effect is subdominant com- 
pared to the first one: the flux pdf is clearly more skewed in 
cDE cosmologies than in the standard case since the cosmic 
web appears to be more evolved in the former case rather 
than in the latter. The shaded area represents the ±la sta- 
tistical error bars inferr ed from a jack-kn ife estimate of high- 
resolution UVES data (|Kim et ahlbOOTJ ) at z = 2.94. 

The flux power is shown in the right panel and we can 
see that a similar (although smaller in magnitude) trend to 
the matter power spectrum case is present. The flux power 
is tilted with a pivot scale roughly corresponding to that of 
the matter power. The differences consist in a suppression 
at the 10% level at the largest scales for RP5 (/3 = 0.2) and 
2% for RP2 (/? = 0.08): the effect is nearly symmetric w.r.t. 
the pivot scale and at 0.01 s/km the power has increased 
by roughly the same amount. The shaded area indicates for 
quantitativ e comparison the ± la statistical error estimated 
at z = 3 bv lMcDonaidI l|2006l ) from the low resolution SDSS 
QSO data set 

If we compare the effects that cDE cosmologies have 
in terms of flux pdf and flux power with those of other pa- 
rameters we can note the following: i) as for flux pdf we 
can co mpare our flndings with Figs. 1 and 2 of lBolton et al.l 
((2003): the effect is similar to that of having a larger value 
for (Tg and a lower value for the parameter 7 for the IGM 
temperature-density relation, but the redshift dependence is 
quite different in the two cases since we have a strong evo- 
lution between 2 = 2 and z = 3 in cDE cosmologies and a 
much smaller one in ACDM; ii) the effect on the flux power 
is distinct from that of having a different as value for the 
30wer spectrum am plitude (see Fig. 13 in McDonald et al.l 



J iower spectrum am plitude (see -big. 13 m McUonald et al.l 
2005h and Fig. 3 in lviel fc Haehneltl (|2006l )): iiJ) the effect 



on the flux power is also different from that of a different 
spectral index which does not allow for a change of sign but 
is either positive or negative in the whole range of scales 
(Fig. 3 in I Viel fc Haehneltl [200^ . To give a rough quanti- 
tative reference for a direct comparison we note that the 
effect in terms of flux power of a non-zero coupling, for the 
particular models investigated here, is similar to that intro- 
duced by a ch ange in the slope of th e linear matter power 
spectrum (see iMcDonald et al.l |2005| ) of about An = 0.05, 
n being the slope of the power spectrum at A; = 0.009 s/km 
and z = 3, while the trends are different for the other cos- 
mological and astrophysical parameters. As for the flux pdf 
the effects are similar to those introduced by a different ag, 
value: for fluxes of the order F=0.9 the effect on the flux 
pdf is around ~ 10% when the erg value is increased by 13%, 
very similar to those introduced by the RP5 model but with 
a different redshift dependence. We also note that the signa- 
ture is degenerate with that of non-gaussianity at the level of 
flux pdf but not at the level of flux power (|Viel et al., ;2009l). 
Therefore the cDE signature on Lyman-a flux statistics for 
the specific models investigated here seems to be unique and 
competitive constraints can be expected. 



3.2 Constraints on the coupling 

We now seek constraints on the coupling parameter /3 by us- 
ing a second order Taylor expansion metho d for t he flux pdf 
and flux power. This meth o d is d escribed in lViel fc Haehneltl 
(2006) and in iLidz et all l|2009l ) and although having the 
drawback of underestimating the error bars it has the ad- 
vantage of being calibrated on accurate full hydrodynamical 
simulations. We perform a Monte Carlo Markov Chain anal- 
ysis in the cosmological and astrophysical parameter space 
by varying all the parameters that impact on the fiux statis- 
tics. We have the following set of cosmological parameters: 
Ho ,ns, as, ^m , a parameter describing the effect of reioniza- 
tion, and the coupling /3; in addition to these we have the fol- 
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lowing astrophysical parameters: T^a (amplitude and slope 
at 2: = 3), 7 (amplitude and two slopes at 2: < 3 and z > ?>, 
To (amplitude and two slopes at z < 3 and z > 3), where 
both the latter functions are parameterized as ^4 [1 + z/A\^ , 
and a parameter describin g the errors indu ced by continuum 
fitting on the flux pd f (see Viel et al.ll2009l ). We use two data 
sets: the flux pdf of iKim et al.1 ((20071) which consists of 63 
data points (21 data points per redshift bin) at z = 2.07, 2.52 
and z = 2.94) and the 132 data points a t z = 2.2, ..., 4.2 of 
the SDSS flux power of lMcDonaldl (|2006l ') (12 measurements 
of the flux power at 0.00141 < k (s/km)< 0.01778 for the 11 
redshift bins) . We include the data covariance matrix in the 
analysis and present results in terms of marginalized values 
for j3. For the joint analysis of pdf and flux power we have a 
total of 15 parameters that are varied and we apply a weak 
prior on TcS (amplitude and slope). There are also other 
13 parameters that describe noise and resolution properties 
and the presence of damped Lyman-a systems for the flux 
power. 

We summarize here the constraints found. For the pdf: 
/3 = 0.08 ± 0.05 and [3 < 0.19 (2cr C.L.) by using the flux 
pdf alone in the range F — [0.1 — 0.8], with a reduced 
X^/f = 1.09 (35 d.o.f.); /3 = 0.04 ± 0.04 and P < 0.1 
(2cr C.L.) by using the flux pdf alone in the whole range 
F = [0 - 1], with a reduced xV^^ = 1-21 (53 d.o.f.). The 
ranges at low and high transmissivity are those that are most 
difficult to model due to the presence of strong systems and 
continuum fitting errors, respectively. Thus, we regard the 
first result presented as more conservative even though we 
do model cont inuum fitting er rors and correct for numeri- 
cal resolution (|Viel et al.ll2009r ). For the flux power only we 
obtain: /3 = 0.07 ± 0.04 and /3 < 0.14 (2ct C.L.) using aU 
the 132 data points {x^/u = 1.16, for 120 d.o.f.). All these 
numbers are reasonable and demonstrate that the regions 
of high transmissivity have a constraining power which is 
stronger than the power spectrum alone. If we c ombine the 
two m easurements we flnd the same trends as in lViel et al.l 
l|2009r ): there is not a very good flt to the data (x^ = 200 
for 164 d.o.f.), and a reasonable x^ is obtained only when 
neglecting the three highest redshift bins of the SDSS flux 
power. In this case, we obtain /3 = 0.05 ± 0.03 and /3 < 0.1 
(2a C.L.) with a reduced x^/i' = 1-09 (146 d.o.f.). AU the 
other parameters are not afl'ected significantly by the new 
parameter introduced and there are not strong degeneracies 
for /3. 

From the analysis performed we can conclude that ro- 
bust 2a upper limits on the coupling constant are in the 
range fi < 0.1—0.2 (depending on the subset of data chosen). 
These bounds are exclusively derived by the analysis of the 
observed properties of the IGM and represent a completely 
new and independen t test of cDE cosmologies w.r.t. previous 
constraints (as e.g. iBean et al.ll2008l : iLa Vacca et al.ll2009l : 
IXial 2009). We regard a 2a limit of /3 < 0.15 as a conser- 
vative overall bound once the statistical limitations of the 
different samples are taken into account. 



4 CONCLUSIONS 

In this work we have explored the possibility of constraining 
the coupling /3 between CDM and DE through the statistical 
properties of the transmitted fiux in Lyman-a forest QSO 



spectra at 2 = 2 — 4.2. For this purpose, we have performed 
the first high-resolution hydrodynamical simulations with 
gas cooling and star formation in the context of cDE models 
and quantitatively exploited the capabilities of fiux 1-pt and 
2-pt functions to constrain the strength of the coupling /? 
between DE and CDM. 

The main results can be summarized as follows: 

• A non-zero coupling between dark matter and dark en- 
ergy produces differeces in the probability distribution func- 
tions for the gas that are more prominent in underdense re- 
gions: the pdf being skewed towards voids in cDE models 
compared to standard cosmologies; 

• Voids in cDE cosmologies are emptier and contain less 
neutral hydrogen as compared to ACDM, and this effect 
might alleviate tensions between simulations and observa- 
tions in voids; 

• The matter power spectra, for the specific cDE models 
investigated here, are also affected in a scale-dependent way 
at high redshift to a level that can be constrained by the 
data; 

• The impact of these effects in terms of flux pdf and 
power spectra is smaller but still larger than present sta- 
tistical errors and more importantly is not degenerate with 
that of other parameters that act in a similar way due to its 
different redshift evolution; 

• By using a Monte Carlo Markov Chain scheme that 
allows to vary all the parameters involved (astrophysical, 
cosmological and noise-related) we obtained a robust and 
conservative upper limit of /3 < 0.15 at the 2a C.L., after 
having marginalized over the other parameters. This limit 
on the coupling /3 is a new and completely independent con- 
straint with respect to previous bounds based on different 
observables. 

This work quantitatively shows that the Lyman-a range 
of scales and redshifts, where the growth of structures can 
be radically different from that measured from a naive ex- 
trapolation of either local or very high redshift probes, is 
promising for constraining coupled dark energy cosmologies. 
The increasing number of QSO spectra that are being col- 
lected (e.g. BOSqj, X-Shootecl) offers the exciting prospect 
of further improving the numbers and of understanding in a 
more refined way, by performing simulations and by address- 
ing systematic errors, the impact that coupled dark energy 
cosmologies can have on the diffuse gas at high redshift. 
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